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We study transport of fermions in a system composed of a short optical lattice connecting two 
finite atomic reservoirs at different filling levels. The average equilibration current through the 
optical lattice, for strong lattice-reservoir coupling and finite temperatures, is calculated within the 
Landauer formalism using a nonequilibrium Green's functions approach. We moreover determine 
quantum and thermal fluctuations in the transport and find significant shot-to-shot deviations from 
the average equilibration current. We show how to control the atomic current by engineering specific 
optical lattice potentials without requiring site-by-site manipulations and suggest the realization of 
a single level model. Based on this model we discuss the blocking effect on the atomic current 
resulting from weak interactions between the fermions. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices have been proven 
to be perfectly suited for implementing physical mod- 
els of interest to the field of atomic and condensed mat- 
ter physics [H, 0| • Specifically, an important part of the 
related experimental efforts have improved and are still 
extending our understanding of nonequilibrium quantum 
transport. These efforts have resulted in the observa- 
tion of fundamentally interesting quantum mechanical 
processes such as Bloch oscillations J3(], Landau-Zener 
tunneling [4J, and interaction-controlled transport |5|. 
Studying nonequilibrium transport in optical lattices has 
several advantages over conventional condensed matter 
systems: Ultracold atoms exhibit slow coherent quan- 
tum dynamics (with kilohcrtz tunneling rates) and are 
detectable in small numbers on microscopic scales. In ad- 
dition, the ability to tune the interactions between atoms 
via Feshbach resonances makes it possible to investigate 
transport of interacting and noninteracting particles. 

While many transport-related experiments employed 
tilted lattice potentials, alternative setups for study- 
ing nonequilibrium phenomena have been suggested, in 
which a current of particles flows between two atomic 
reservoirs. Micheli et al. |fjH8| pointed out that, analo- 
gous to a transistor, an impurity atom may be utilized 
to control the flow of a one-dimensional Bose or Fermi 
gas. They considered the full time-dependent coherent 
dynamics of the ultracold gas in a closed system and de- 
termined the particle current by using analytical approx- 
imations [7j and time-dependent density matrix renor- 
malization group (DMRG) calculations Q- 

Pepino et al. (jj, [l(| generalized this idea by replac- 
ing the single impurity with an optical lattice coupled 
to separate bosonic reservoirs in order to emulate the 
behavior of semiconductor electronic circuits (see also 
[ill]). In Refs. [9|, [To] reservoirs were introduced and 
specified as large sources or sinks of particles at zero 
temperature with a Fermi-sea-like energetic distribution, 
constant chemical potential, and fast-decaying system- 
reservoir correlations. A quantum master equation, re- 
lying on weak system-reservoir coupling, was used to de- 



scribe time evolution of the system, thereby in part elim- 
inating the coherent evolution of the reservoir. 

Here we consider the evolution of a one-dimensional 
Fermi gas loaded into an optical lattice which is parti- 
tioned into two large incoherent reservoirs L and R con- 
nected by a short coherent region C, as illustrated in 
Fig.[TJ In this setup the difference in the chemical poten- 
tials fir, and hr of reservoirs L and R drives a current of 
fermionic atoms through the coherent region C . 
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FIG. 1. (Color online) Fermions confined to an optical lat- 
tice hop from the left reservoir L through a short coherent 
part C (sites 1 to m) into the right reservoir R. The hopping 
parameters Jij, the on-site energies Ej, and the couplings Jl 
and Jr may take arbitrarily engineered values, whereas the 
hopping parameter in the reservoirs Jo is held constant. 



In contrast to previous theoretical works, we take a 
mesoscopic perspective on nonequilibrium transport be- 
tween the atomic reservoirs: We analyze the evolution 
of the Fermi gas within the Landauer formalism [12L 1 131 ] , 
where transport is described as a transmission process 
through the coherent region C with fermions emitted and 
absorbed by the reservoirs L and R, respectively. In this 
vein, we understand reservoirs to be finite containments 
filled with noninteracting fermions in the ground state, 
characterized by a Fermi distribution with finite tem- 
perature and well-defined (but not necessarily constant) 
chemical potential. The transmission is determined by 
a nonequilibrium Green's functions (NEGF) approach in 
the tight-binding picture [Tj, [l5| . This approach allows 
us to consider a wide range of experimentally relevant 
parameters. Moreover, we express our results in terms of 
the number of atoms accumulated in the reservoirs, which 
(in particular for small currents) may be more accessible 
in experiments than the current. 



Accordingly, we com plem ent and significantly extend 
the results in Refs. [9|, [l0| in several directions for the 
case of fcrmionic atoms: First, we take into account that 
experimental ultracold systems arc finite and therefore 
we describe not only instantaneous steady-state currents 
at constant chemical potentials but the full equilibration 
process between the reservoirs. Second, our analysis is 
valid for arbitrarily engineered configurations of the co- 
herent region C and for strong coupling to the reservoirs, 
thereby going beyond the weak-coupling expansion in 
Refs. [a, G3- Strong coupling indeed corresponds to the 
most elementary setup with the parts L, R and C con- 
sisting of a single homogeneous optical lattice. Third, by 
using the tools of full counting statistics [Tj] and consid- 
ering intrinsic damping mechanisms we explicitly evalu- 
ate the evolution of the fluctuations in the current. These 
fluctuations indicate large shot-to-shot deviations from 
the average value in repeated measurements of the num- 
ber of particles in the reservoirs. Our approach includes 
finite-temperature effects relevant to both the average 
current and the fluctuations. For instance, thermal fluc- 
tuations are shown to build up during the equilibration 
process until they reach a constant value proportional to 
the temperature of the Fermi gas. 

In the first part of this paper we analyze the evolution 
of the Fermi gas within the NEGF-Landauer formalism. 
In this framework we discuss the properties of the system 
in terms of the average current, the filling levels of the 
reservoirs, and fluctuations of these quantities. In the 
second part we first apply the formalism to the case of 
constant transmission between the reservoirs and subse- 
quently discuss more complex situations. In particular, 
we show how to control the current by modulating the 
connecting optical lattice with the help of additional op- 
tical potentials. Finally, we suggest a way to realize a 
single-level model and demonstrate that in this model 
even weak interactions between the fermions are suffi- 
cient to reduce the atomic current. 



and Jo, respectively, and Ej are on-site energies. The cou- 
plings Jl and J R connect the central part to the reser- 
voirs which are each composed of M sites. 

The Hamiltonian Hs, not specified explicitly, repre- 
sents interactions of the reservoirs with an engineered 
environment, e.g., with an atomic gas or optical radi- 
ation that is not necessarily far detuned. The interac- 
tions introduced by Hs add dissipative and incoherent 
processes to the reservoirs so that they act as scmiclassi- 
cal systems equivalent to metallic electrodes. These pro- 
cesses are assumed to destroy coherence and to relax the 
fermions to the ground state on a time scale shorter than 
h/ Jo, making it possible to attribute a Fermi distribution 
with well-defined temperature and chemical potential to 
the reservoirs. An explicit scheme suggested in Ref. [16[ 
achieves this aim through a combination of coherent laser 
excitations and dissipation into an ambient supcrfluid. 
First, fermions with high momentum arc transferred into 
the first excited band of the optical lattice via a Raman 
process. Subsequently, the excited states decay into the 
lowest Bloch band due to emission of phonons into the 
supcrfluid. An iteration of this procedure results in a 
stable Fermi distribution of the atoms in the reservoirs. 

To apply the Landauer formalism modified to account 
for the finite size of the reservoirs, we introduce the num- 
ber operators N a , with a = L or R. measuring the num- 
ber of fermions in the reservoirs and the expectation val- 
ues (N a ). We specify the state of the system by the aver- 
age particle number (N a ) and the current dt{N a ) through 
the central part. To obtain the current in the Landauer 
formalism we treat the connecting optical lattice as a 
scattering potential with the energy-dependent transmis- 
sion T(e). Hence the average current dt{N R ) = —di{Ni) 
is the sum of all possible scattering transfers between the 
two reservoirs 



dt(N R (t)) 



de 
2nh 



T(e)[f L (s,t)-f R (e,t)} 



(2) 



II. NEGF-LANDAUER MODEL 

We start with the theoretical framework required to 
determine the average equilibration current through the 
lattice and the filling levels of the reservoirs. The Hamil- 
tonian of the system within the lowest Bloch band and in 
the tight-binding approximation is He + Hl + H R + Hi 
with 



H c = - Yl J m4cj + Yl £jc]cj h 3 e l.-m 

(id) 3 

H L =H R = -JoY J eje,- + H S i,j<£l..m (1) 

(id) 

Hi = -Jl(c\co + cjci) - JR(c^c m+ i + c\ n+1 c m ) , 

where the central part C is formed by the sites 1 to m and 
(i,j) denotes the sum over nearest neighbors. The oper- 
ators c*- (£j) create (annihilate) a spin-polarized fermion 
in a Wannier state localized at site j. The hopping pa- 
rameters in the central part and in the reservoirs are Jjj 



The Fermi functions of the reservoirs are f a (e,t) = 
[ e (e-M Q )/fcsT _|_ i] — i^ w ith the Boltzmann constant k B , 
the temperature T, and the time-dependent chemical po- 
tential n a (f). Note that Eq. ([2]) is approximately valid 
provided that f a (s, t), or equivalently jx a {t), varies slowly 
on the microscopic time scale H/Jo- This can readily be 
achieved by either increasing the size of the reservoirs 
M or decreasing the transmission through the central 
part C. 

The chemical potential \x a is related to the particle 
number (N a ) since the reservoirs are finite and the total 
number of fermions in the system is fixed. The implicit 
relation between [i a and the particle number (N a ) is 



(N a ) = M dep(e)f a (e), 



(3) 



with the density of states of the reservoir given by p(e) = 
1/tt^(2J ) 2 - (e-2J ) 2 for < e < 4J and zero other- 
wise. Upon solving Eq. ^ for fj, a (t), either analytically 
or numerically, one finds the chemical potential \i a (t) as a 
function of (N a (i)). Thus, as a consequence of the finite- 



size reservoirs the Landauer formula in Eq. ([2]) becomes 
a closed (intcgro-diffcrential) equation for (N a (t)). 

The transmission T(e) is efficiently determined by use 
of the NEGF approach [Tj, [l5( . In terms of Green's func- 
tions we have 



T(s) = \G hm (e)\ 2 T L (e)r R (£), 



(4) 



where Gij(e) is the full retarded Green's function for 
the central part and T a (e) = 2irJ 2 p a (e) describes the 
coupling to the reservoirs. Here, pl(e) and p R {e) arc, 
respectively, the local density of states at sites and to+1 
(see Fig. QJ. Consequently T{e) encodes the coherent 
evolution of the central part, governed by He, as well as 
the coupling to the reservoirs. 

To calculate the full Green's function Gij(e) from He 
we start with the bare Green's function gij(s) whose 
components obey the equation 



(e - ei)9i,j + V Ji,k9k,j = 5; 



k 



'i -J 



or equivalently g(e) = l/(e — He), with both the Green's 
function and He in matrix notation. The bare Green's 
function gi,j(e) with poles at the energy levels of He 
would be sufficient to determine the transmission if the 
coherent region C were coupled weakly to the reservoirs. 
However, the full propagation of a fermion between the 
sites i and j includes excursions into the reservoirs due to 
the coupling, which leads to the broadening T a ~ J 2 / Jo 
of the energy levels. More precisely, the excursions re- 
sult in corrections to <?ij(e) in terms of self-energies 
s £ = j l9o,o and Y, R = J R g m +i,m+i, where <fc,j(e) de- 
note the Green's functions of the reservoirs. Including 
the coupling to the reservoirs to all orders we arrive at 
the Dyson equation for the full Green's function 



Gi,j — gij + ffi,lSiGl,j + gi,mJ-iRGm,j 



(5) 



The relevant matrix element for the transmission G\. m 
is then obtained from Eq. ([5]) by solving a set of simul- 
taneous equations for Gi.i, G m . m , G-\ m and G m .\- As a 
result one finds G\, m = g\. m jD with [14| 

D = (l- S i gi i i)(l - £ fl g m , m ) - Y, L Y, R gl m 

assuming that gi Jn = g m ,i- 

For our specific setup we treat the reservoirs as semi- 
infinite optical lattices, for which the Green's function 
(jij at the end site reads [15| 



ffo.oOO = [(e - 2J ) - V(2J ) 2 - (e - 2J ) 2 J /2 J 2 

and 5m+i, m +i = <7o,o; the local density of states is given 
by p a {e) = — (l/7r)Imgo,o(£)- These relations allow us 
to find explicit expressions for the couplings r Q and the 
self-energies E a in Eq. ((5). 



III. FLUCTUATIONS AND DAMPING 

So far our analysis was restricted to the average 
fermion number (N a (t)) and the current dt(N a (t)), 
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FIG. 2. (Color online) Filling level and fluctuations in 
the right reservoir (initially empty) for a constant transmis- 
sion 7o. The average filling level {Nn{t))/M (solid line) in- 
creases with time t/t eq and saturates at the equilibrium value 
1/2. The standard deviation y/{5N^(t))/M from the average, 
due to thermal fluctuations, is indicated by the orange (gray) 
band. The Fano factor F = (SN%(t))/{N R (t)) (dashed line) 
decreases with time and approaches a constant value in the 
regime t/t eq 2> 1. The parameters are /csT/JoM = 1/10 for 
the filling level and IcbT/Jo = 1 for the Fano factor. 



which are found from Eqs. © and ([3]) for a given trans- 
mission T(e) in Eq. Q. We now turn our attention to 
quantum and thermal fluctuations present in the system. 
To simplify the problem of determining the fluctuations 
we treat creation and damping processes separately and, 
if possible, add their effects together. 

Fluctuations created during the evolution of the sys - 
tcm can be found by using the Levitov formula (l3j |. 
which yields the full counting statistics of thermal and 
quantum fluctuations provided (N a (t)) and hence fj. a (t) 
are known. However, we limit our analysis to the most 
relevant statistical parameter, namely the variance of the 
number of fcrmions (5N%) = ((N a - (N a )) 2 ). With ini- 
tially no fluctuations present we find for the variance 
(8N 2 (t)) after a time t according to the Levitov formula 

(5N*(t)) =j\ s j^L {T(e)t(e)[f L (e, s) - f R (e, s)} 2 



T(e)[f L (e, S )f L (e, S ) + f R (e, S )f R (e,s)}} 



(6) 



where we introduced f a (e,t) = 1 — f a (s,t) and T(e) = 
1 — T(e). At zero temperature, Eq. © describes the cre- 
ation of quantum fluctuations caused by the probabilistic 
nature of the particle transfer through the optical lattice. 
On the other hand, at equilibrium between the reservoirs 
and for T(e) = 1 the fluctuations are purely thermal. 

Unlike in conventional mesoscopic systems with 
infinite-sized electrodes, the fluctuations described by 
Eq. ([B]) are constantly damped out at a rate 7. This 
intrinsic damping occurs because fluctuations in the cur- 
rent immediately lead to fluctuations of the chemical 
potentials, which drive the system back to the evolu- 



tion according to the mean- field description (N a (i)). 
To find an explicit expression for the intrinsic damp- 
ing 7 we determine the change in the chemical poten- 
tial 5 n a caused by an excess of particles 6N a -C (N a ) 
with respect to (N a ). From Eq. ^ we obtain (N a ) + 
SN a = M J dep(e)f a (e,p a + Sp, a ), which to lowest or- 
der in 8p a and SN a results in the linear dependence 
Sfi a = SN a /MpT([J-a) with 



Pt(Po 



i^r/ d ^ (£)scch2 



2k B T 



(7) 



Note that pT{p a ) reduces to p(p a ) in the case T = 
and to k R T /A in the limit of infinite temperature. Sim- 
ilarly, we expand the Landaucr formula to lowest or- 
der in 5p a and dt(SN a ) and use the relation 5p a — 
8N a /MpT{p a ) to find the time dependence for small fluc- 
tuations dt(SN a ) = —"f5N a with the damping factor 



1 



7 = 



2-irhM 



Tt(pl) . T T (p R ) 



Pt(pl) Pt(pr) 



(8) 



where Tr(Pa) is defined in the same way as pT(p a ) in 
Eq. ([7]) . The factor 7 depends on the filling level of the 
reservoirs and is therefore time dependent in general. Im- 
portantly, 7 is always positive and thus fluctuations are 
indeed damped out. 

As a result, the fluctuations in the system depend not 
only on the properties of the coherent region but also on 
the reservoirs via their density of states p(e) evaluated at 
the filling level. Thus by choosing the appropriate reser- 
voirs it should be possible to perform experiments in ei- 
ther the fluctuation-dominated or the mean-field regime. 
This works particularly well for an optical lattice reser- 
voir since its density of states varies considerably over 
the entire bandwidth 4 Jo so that the intrinsic damping 7 
can be tuned over a wide range. 



IV. CONSTANT TRANSMISSION 
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FIG. 3. (Color online) The engineered transmission T(e) as 
a function of the energy e/Jo for different modulations of an 
optical lattice with length m = 10. (a) A single centered beam 
with waist a — 2 and increasing depths V/ Jo = 1, 2, 4 (dotted, 
dashed and full line) shifts the energies out of the reservoir 
band and reduces the transmission to a peak at the upper 
band limit, (b) Two beams isolate a few central lattice sites 
and create a single resonant level at e/Jo ~ 1 for sufficiently 
strong intensities. The beams are positioned at v\ — 3, ^2 = 8 
with waist a = 1, depths V2 = 1, 3, 5 (dotted, dashed and full 
line) and V1/V2 = 1/2. 



To gain physical insight into the transport between 
reservoirs we apply our general results to the impor- 
tant special case of constant transmission: T(e) = 7o 
for < e < 4Jo and zero otherwise, where the con- 
stant transmission 7o takes values < 7o < 1- In this 
case the Landaucr formula for the current reduces to 
dt(Nfi) = ToAp/2irh and hence one only has to deter- 
mine the dependence of Ap = p,i — pr on the average 
particle number (N a ) and the temperature T. 

For concreteness we consider the equilibration process 
in the low-temperature regime fcsT <C Jo, neglecting cor- 
rections of the order (fc^T/ Jo) 2 . In this case the chemical 
potential difference is Ap = 4 J cos(ir {Nr)/M) assuming 
the initial conditions (Nl) = M and (Nr) = 0. Solving 
the Landauer equation we obtain the evolution of the 
filling level 

(N R (t)) = CU &rct&n[smh(t/t eq )] 

and the current through the central region dt(Nu(t)) = 
(U/R) sech(t/t eq ). The equilibration time scale t eq = RC 



satisfies t eq ^> H/Jq in the relevant parameter regime 
as required for Eq. ([2]) to be valid. Here, analogous 
to a classical RC circuit, we introduced the resistance 
R = 2irh/To, the capacitance C = M/AttJq, and initial 
capacitor voltage U = 4J . 

We next determine the evolution of the fluctuations 
during the equilibration process. For simplicity we as- 
sume from now on that 7o is close to unity, i.e., 7o ~ 1, so 
that predominantly thermal fluctuations with (6N%(t)) = 
2tksT/R are created during the equilibration process, 
as can be found from the Levitov formula in Eq. ^. 
Since correlations of thermal fluctuations decay fast on 
the time scale t eq , the evolution of the deviations SNr 
from (Nn(t)} can be expressed in the form of a Langevin 
equation 

d t {5N R ) = - 7 (t)SN R + t;(t). 

The time-dependent damping factor resulting from 
Eq. © reads 7 (t) = sin[n {N R (t))/M]/RC and thus 



damping increases as the system approaches equilib- 
rium. The fluctuations are represented by the sto- 
chastic force £(£) and satisfy the condition (£(£)£(£')) = 
(2ksT/R) 5(t — t '). Using standard techniques for sto- 
chastic problems [17| we nn d the evolution of the fluctu- 
ations 

(SN R (t)) = Ck B T [scch 2 (t/t eq ) t/teq + taok{t/teq)] ■ 

We see that the fluctuations increase linearly in the 
regime t/t eq <C 1, where damping is weak according to 
Eq. (|SJ). On the other hand, in the limit t/t eq ~^> 1 the 
fluctuations converge to the constant value CksT, which 
results from the competition between thermal fluctua- 
tions and the intrinsic damping 7 = 1/RC at equilib- 
rium. 

Figure [21 shows the time evolution of the average filling 
level (Nn(t))/M and the fluctuations in the right reser- 
voir. The filling level increases linearly in the regime 
t/teq <C 1 and saturates at the equilibrium value (Nr) = 
(Nl) = M/2. Thermal fluctuations around the aver- 
age are limited by damping; however, they indicate sig- 
nificant shot-to-shot deviations from the average filling 
level. Of particular interest is the Fano factor defined 
by F = (5N R (t))/(N R (t)), which is independent of the 
size of the reservoirs M and gives a direct measure of the 
temperature of the system. The Fano factor converges to 
the value F = kBT/2irJ as the system equilibrates, i.e., 
in the limit t/teq ^> 1- This result even holds for arbi- 
trary constant transmissions 7o as quantum fluctuations 
are damped out in this limit. 

An elementary experimental configuration with con- 
stant transmission, namely 7o ~ 1, consists of a sin- 
gle homogeneous optical lattice partitioned into a suf- 
ficiently long coherent part C and the reservoirs L and 
R. The corresponding hopping parameters and couplings 



0.4 



are J, 



Jo and J a = Jq, respectively. In fact, the 



strong couplings J a lead to significant broadening of the 
cosine-distributed energy levels of the coherent part [18| ■ 
As a consequence, the broadened energy levels merge to- 
gether in the regime J a 3> 27rJo/m, which results in an 
approximately constant transmission. 
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FIG. 4. (Color online) Filling level and fluctuations in the 
right reservoir (initially empty) for a lattice modulated by a 
single beam, (a) The average filling level (Nn(t))/M (solid 
line) as a function time t/t eq (with t eq = HM/2Jo) for the 
modulation strength V j Jo — 2. Different bands indicate the 
standard deviation sj (5 N^(t)} j M due to quantum fluctua- 
tions (blue or light gray) and total fluctuations, i.e., quantum 
plus thermal (orange or dark gray), (b) The same quanti- 
ties for a stronger modulation strength Vj Jo = 4. The lower 
transmission results in a smaller current between the reser- 
voirs and enhanced quantum fluctuations in comparison to 
(a). In both plots damping has not been taken into account 
and the parameters are fcsT/Jo = 1, m — 10 and M = 50. 



V. TRANSMISSION ENGINEERING 

The usefulness of the NEGF-Landauer formalism is 
most evident if we want to calculate the current through 
an engineered optical lattice with tailored parameters Jij 
and £j. Since the optical lattice potentials for the reser- 
voirs and the central region would most likely be pro- 
duced by the same counterpropagating laser beams we 
set Jij = J a = Jq and focus on modified on-site energies 
By A possible experimental configuration with tailored potential shifts the energies Ej out of the reservoir band 



with the potential strength V and width a measured in 
units of lattice spacings. Depending on the detuning, V 
may take positive or negative values 19| . We emphasize 
that unlike the scheme proposed in Ref. 9| such a config- 
uration does not require site-by-site control of the optical 
lattice, neither of the on-site energies Ej nor the hopping 
parameters J,j or J a . 

A specific setup consists of a single laser beam, cen- 
tered at v w m/2 and with beam waist a ~ m/A. acting 
as a V^-dependent switch for the fermion current. The 



on-site energies Ej involves one (or several) laser beams 
crossing the central region perpendicular to the optical 
lattice. The optical potential caused by a single beam 
centered at position v shifts the energies as 

Ej = Vcxp[-(j -vf/a 2 ], 



and thus reduces the transmission significantly, as shown 
in Fig. G3a). This configuration provides the possibility 
to study the dependence of the fluctuations on the trans- 
mission T"(e) and the temperature T of the fermions: If 
we choose the initial conditions (Nl) = M, (Nr) = 
and stay far from equilibrium then according to Eq. (|5J) 



damping is negligible. As a consequence, thermal and 
quantum fluctuations lead to significant deviations of the 
filling level from (N R (t)), which are detectable by count- 
ing the actual number of fcrmions in the right reservoir. 
Figure [4] shows the average particle number (N R (i)) and 
the expected fluctuations for modulations of the optical 
lattice with two different potential strengths V. At zero 
temperature, only quantum fluctuations caused by the 
limited transmission contribute, whereas at finite tem- 
perature fluctuations are further increased. By compar- 
ing Figs.Hta) and|4|b) we see that quantum fluctuations 
become more important for reduced transmissions, i.e., 
for stronger potential strengths V. 

A paradigmatic system in the context of mesoscopic 
physics is the single- level model [la], or in the case of 
interacting fermions the Anderson impurity model [201 ] . 
These models may be used, e.g., to study the Kondo ef- 
fect or to describe transport through a single quantum 
dot. Realization of a single- level model can be achieved 
by means of two laser beams with different detuning lead- 
ing to energy shifts 

e, = V 1 exp[-(j - vi?l° 2 \ - V 2 exp[-{j - v 2 ) 2 /a 2 ] . 

If the beams are separated with 1^1 — ^2 ~ to/2 and 
narrow a ~ m/4 then for sufficiently strong potentials 
Vi , V2 > the transmission exhibits a single peak, as 
shown in Fig. EJb). The position of the peak depends on 
the ratio V1/V2 and the strength of the potentials deter- 
mines the width of the peak, i.e., the effective coupling 
to the reservoirs. 

The emergence of the single level can be understood 
in the energy band picture: The first beam shifts the 
unmodulated band of width 4 Jo upward, while the sec- 
ond beam shifts the band downward. As a consequence, 
a small region between the beams is isolated from the 
reservoirs and acts as a single energy level. The effective 
couplings T Q of the level to the reservoirs, or equivalcntly 
the width of the single level, is readily controllable by the 
strengths of the beams V\ and V 2 ■ This makes it possible 
to access the weak-coupling regime considered in Ref. @ 
as well as the strong-coupling regime without the require- 
ment of specific control of the hopping parameters J a . 



VI. EFFECT OF INTERACTIONS 

Let us now discuss the effect of weak interactions be- 
tween the fermions on the basis of the single-level model. 
If we consider a spin-balanced mixture of fermions in 
two different internal states, denoted by up and down, 
then the total Hamiltonian in Eq. ([1) is augmented by 
the interaction term Hjj = cfVnJnj. Here, U is the 
either positive or negative interaction strength between 
the fermions and nj (n|) is the occupation number op- 
erator for the up (down) states. For our analysis we use 
the mean-field approximation fv.h, « nj(n^) + (fij)h^ — 
{fij){n^) which we expect to be valid for sufficiently small 
reservoirs and low filling levels [2l| . The chemical poten- 
tials of the reservoirs /i a are changed accordingly due to 
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FIG. 5. (Color online) The effect of interactions on a single 
level with original (U = 0) position £0 = and the relative 
chemical potentials (il/Jo = 1 and [ir/Jo = — 1. The occu- 
pation (n) (blue or dark gray) and the energy shift U{n)/Jo 
(orange or light gray) as a function of the interaction U/Jo 
are plotted for the couplings T/Jo = 0.1 (solid) and T/J = 1 
(dashed). For U > the level is depleted and shifted toward 
jiL- For U < the level is almost completely occupied and 
shifted below [j,r for sufficiently strong interactions. The cur- 
rent is blocked if the single level at eo + U{n) leaves the energy 
window [(j,l , A*i?] ■ 



the interactions; e.g., the chemical potential for spin- up 
states is given by /xj, = fx a H- U (n^). In the following we 
will however neglect these trivial changes. 

The effect of interactions on the single level is to shift 
its original (U = 0) position £0 by the interaction energy; 
e.g., for spin-up fermions the level position is £0 -I- U(h^). 
The corresponding Green's function reads 

&(e) = l/[e - £ - U(&) + T L (e) + T R (e)} (9) 

with the effective couplings T a . The occupation of the 
single level is given by 



(10) 



<^> = / % |G%)| 2 [r L (£) f L ( £ ) + T R (e)f R (e)] , 

with the time dependence of the Fermi functions omitted. 
In the case of a spin-balanced mixture with (n^) = (n^) 
the Green's function C~(e) and the occupation (n^) can 
be determined sclf-consistcntly to obtain the transmis- 
sion T(e) oc |G^(£)| 2 , which is the same for both internal 
states. In principle, the average particle number (N a (t)) 
and the current dt(N a (t)) are then evaluated as for non- 
interacting fermions. 

The effect of interactions on the transmission can be 
qualitatively understood in the wide-band limit with con- 
stant couplings r Q = r. Inserting G^(e) into Eq. (fTU)) 
then yields the self-consistent equation for the occupa- 
tion [2CJ 

/-\ 1 , 2 V^ *. ( V* -e -U(h)\ , . 

< n > = 2 + n ^ n I 2T J ' (U) 



where the average occupation number applies to both 
spin states; i.e., (h) = (fV) = (rr - ). Figure [5] shows the 



average occupation (h) of the single level and the energy 
shift U(n) as a function of the interaction strength U ac- 
cording to Eq. (JTTJ . For repulsive interactions U > we 
observe a depletion of the single level and a shift to higher 
energies. This shift is bounded by /i£ in the limit of van- 
ishing coupling r, but takes values larger than /j,i in the 
case of finite T. For attractive interactions U < the oc- 
cupation (h) increases and makes an abrupt transition to 
(n) w 1 accompanied by a shift of the level below /j,r. In 
both cases the single level at e + U(n) eventually leaves 
the energy window between [i^ and \i^ and hence the 
current through the level is strongly suppressed. Thus 
interactions offer an alternative approach to control the 
current through the lattice. 



VII. CONCLUSIONS 

Using the mesoscopic NEGF-Landauer approach we 
have analyzed nonequilibrium transport of fermions 
through an engineered optical lattice for arbitrarily 
strong coupling to two reservoirs at finite temperatures. 
We have characterized the full equilibration process by 
calculating the accumulated number of atoms in the fi- 
nite reservoirs, which is a directly accessible quantity in 
experiments. Considering experimentally relevant sys- 
tem parameters we found that the reservoirs equilibrate 
on time scales comparable to the duration of typical ul- 
tracold atom experiments. Our systematic analysis of 
created and damped fluctuations in the finite system re- 
vealed that the mean-field description gives an incom- 
plete picture of fcrmion transport since significant shot- 
to-shot variations from the average current, partly due 
to thermal effects, are to be expected. This is directly 
relevant to the emulation of semiconductor electronic cir- 
cuits, where preferably single-shot measurements are re- 
quired to determine the current [lOj. 

The study of fluctuations around the average current 
revealed additional information about the processes tak- 
ing place in the system: We found that thermal fluc- 
tuations build up on the time scale of the equilibration 
process until they reach a constant value proportional to 
the temperature of the Fermi gas. As an aside, we note 
that thermal fluctuation between equilibrated reservoirs 



may therefore be used for thermometry of the system. 
On a more fundamental level, we saw that a decrease 
in the current due to a lower transmission of the coher- 
ent region is necessarily accompanied by higher quantum 
fluctuations. This correlation allows the experimenter, 
e.g., to distinguish between changes in either the chemical 
potentials or the transmission as the cause of a reduced 
current. 

We have shown that modulations of a homogeneous 
lattice potential can be used not only to reduce the equi- 
libration current, but also to realize a single- level model 
with full control over the position and the coupling of 
the level. This setup requires neither additional impu- 
rity atoms nor site-by-site manipulations of the optical 
lattice. Moreover, advanced experimental techniques for 
producing tailor-made optical potentials, by employing 
either acousto-optical deflectors [22[ or holographic mask 
techniques [23[ , are expected to further facilitate the cre- 
ation of engineered optical lattices. Finally, our mean- 
field analysis of interaction effects revealed that even 
weak interactions between the fermions suppress the cur- 
rent through a single level, which can be exploited to 
control the current. 

We conclude with the observation that our approach to 
nonequilibrium transport between finite reservoirs may 
be applied to similar ultracold atomic setups j24| or, more 
generally, to mesoscopic systems such as electrons on liq- 
uid helium [25| . Possible extensions of this work include 
transport of bosonic atoms, similar to the analysis in 
Refs. |7J, 126J, quantum pumping between reservoirs using 
time-dependent modulations of the optical lattice [27j |. 
and the effect of interactions between fermions on quan- 
tum fluctuations S3 ■ 
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